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Abstract The appearance of eruptive space plasma processes, e.g., in eruptive 
flares as observed in the solar atmosphere, is usually assumed to be caused by 
magnetic reconnection. The process of magnetic reconnection is often connected 
with singular points of the magnetic field. We therefore analyse the system of 
stationary resistive/non- ideal magnetohydrodynamics (MHD) in the vicinity of 
singular points of flow and flcld to determine the boundary between reconnec- 
tion solutions and non-reconnectivc solutions. We flnd conditions to enable the 
plasma to cross the magnetic scparatrices also inside the current sheet, close 
to the current maximum. The results provide us with the topological and ge- 
ometrical skeleton of the resistive MHD flclds. We therefore have to perform 
a local analysis of almost all non-ideal MHD solutions without a specific non- 
idealness. We use Taylor expansions of the magnetic field, the velocity field and 
all other physical quantities, including the non-idealness, and with the method of 
a comparison of the coefficients, the non-linear resistive MHD system is solved 
analytically. In the vicinity of a stagnation point, it is reasonable to assume 
that the density is constant. Wc find that the electric field has to be zero and 
that the non-ideal term/resistivity has to depend on the spatial coordinates and 
cannot be constant, otherwise it has to be zero everywhere. It turns out that 
not every non-ideal flow is a 'reconnective' flow and that pure resistive/non-ideal 
MHD only allows for such reconnection-like solutions, even if the non-idealness 
is localized to the region around the magnetic null point. It is necessary that the 
flow close to the magnetic X-point is also of X-point type to guarantee positive 
dissipation of energy and annihilation of magnetic flux. If the non-idealness has 
only a one-dimensional, sheet-like structure, only one separatrix line can be 
crossed by the plasma flow, similar to reconnective annihilation solutions. 
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1. Introduction 



Magnetic reconnection is thought to be a process, being responsible for many 
eruptive plasma phenomena in space plasmas and astrophysical plasmas, like 
geomagnetical substorms or eruptive flares. Although magnetic reconnection 
in two dimensions (2D) is fairly well understood, e.g., see the comments in 
|Baty, Forbes, and Priest| (|2009p , it would be interesting to have more detailed 
informations about the topological and geometrical structure of flow and field 
lines in the vicinity of the singular points of plasma flow and magnetic field. 

Classical reconnection scenarios, see iPctschek ( 1964, ) or Sweet-Parker model 
(see, e.g.. lSweet[[T958|) propose a magnetic null point and a stagnation point flow 
into the diffusion region, i.e. the stagnation point is inside this diffusion region. 

[Priest and Cowley] (|1975p where the first to analyse the case of incompressible 
2D MHD with constant resistivity. It turned out that either the magnetic field 
must be of higher order in the spatial variables x and y or that there is no 
stagnation point flow, but a shear flow. Therefore their result is, that the clas- 
sical 'hyperbolic' stagnation point flow needs higher order terms, concerning 
the spatial variables ix and y). Solutions locally containing only higher order 
terms do not allow for topologically/structurally stable magnetic fields. There- 
fore either no reconnection can take place, as in the case of higher order terms 
no 'hyperbolic' magnetic field can exist, or the stagnation point fiow is not of 
hyperbolic type. 

Annihilation solutions have been studied, where [Craig and Henton] p995p 
chose a special ansatz for the solution of the resistive MHD to get reconnec- 
tion solutions. They emphasize the first order momentum equation and neglect 
the energy transport (equation) or rather the entropy conservation (equation), 
starting with a nonlinear perturbation of magnetic annihilation solutions. This 
lead them to so called 'reconnective annihilation' solutions, where only one of 
the two separatrix-lines are crossed, and the other is only tangent to the con- 
verging streamlines. The current sheet has a one-dimensional structure (straight 
line). The results found by [Craig and Henton and Craig and Rickard] (|1994p 
confirm the results found earlier by Priest and Cowley (|1975p , who found more 
'shear-like' flows instead of typical stagnation point flows. 

Later on [Tassi, Titov, and Hornig (2002^ and [Titov, Tassi, and Hornig[ 
P004p extended the method to curvelinear current sheets. 

It was shown by IPriest et al.\ p994p and later on in extended form by 



Watson and Craig (|1998p that under certain circumstances (like constant re- 



sistivity or current depending/anomalous resistivity and sub-Alfvenic flow etc.) 
reconnection is impossible, the so called anti-reconnection theorems. 

In 3D for constant resistivity a careful analysis of topologically different solu- 
tions has been discussed in [Titov and Hornig] (|2000p . In 2D such an analysis is 
missing and should be done here, but without taking the restriction to constant 
resistivity into account. 

In contrast to the aforementioned models we do not search for special solu- 
tions, but the most general solution without constant or special non-constant 
resistivity or a specific non-idcalncss. Thus our analysis covers all forms of non- 
ideal terms. In the sense of topological fiuid dynamics we want to get insight 
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how the field and streamfines are rooted in the null point(s) of flow and field, i.e. 
which geometrical shapes of field and streamlines are possible and which not. 
We define and investigate the infiuence of an 'effective' non-constant resistivity 
and use hereby the full energy equation of resistive MHD instead of using only 
the assumption of incompressibility. 

The problems concerning (exact) analytical models very often are: 

• The physical quantity 'resistivity' is often only a constant smallness param- 
eters 

• if really gradients of the resistivity r/ and maybe also the amplitude/absolute 
value are recognized as important for the non-ideal process, what are the 
'shapes' of such resistivities enabling magnetic reconnection? 

• it is not clear which flow topologies are allowed to generate reconnection, 
reconnective annihilation solutions, or other solutions for general non- ideal 
or resistive terms 

Thus our aims are: 

• To show not every non-ideal or resistive process in the vicinity of a null 
point is a magnetic reconnection process 

• it is not enough to have a localized resistivity to get a reconnective solution, 
finding parameters that mark the boundary between reconnection and non- 
reconnection solutions 

• why should there be no reconnection process close to the magnetic null point 
and what happens if we do not restrict ourselves to complete incompressible 
dynamics without energy transport? 

• finding analytical and exact solutions close to the null point 

• detailed investigation of topological skeleton of flow and magnetic field in 
the frame of resistive MHD; it is an analysis of the resistive system close to 
singular points of flow and magnetic field 

• performing linearization, necessary to get information about the skeleton 
of magnetic reconnection 

We will concentrate on the topology and geometrical properties of fiow and 
field and discuss different cases with respect to their implications for magnetic 
reconnection. 

To summarize: The main goal is to get the topological and geometrical skele- 
ton of the resistive MHD equations, i.e., if the Jacobian of the magnetic field 
at the magnetic null point has only real Eigenvalues (hyperbolic null point), 
which Eigenvalues can be found for the Jacobian of the plasma velocity at the 
stagnation point? In our analysis the null point of the velocity field (=stagnation 
point), should have an almost identical position and a vanishingly small offset to 
the magnetic null point. The resistvity 77 should be a positive quadratic form in 
X and y close to the magnetic null point, indicating that the resistivity indicates 
real dissipation in the form of Ohmic heating and annihilation of magnetic flux, 
see also the explanations and calculations in the section [3?3l How does the geo- 
metrical structure of the velocity field look like, and which kind of Eigenvalues 
are allowed by the resistive or in general non-ideal MHD equations? 
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2. Assumptions and basic equations 



2.1. The topological and geometrical structure of the magnetic field 

The topological classification of 2D vector fields in the vicinity of their null points 
is described in the textbooks (with connection about phase portraits of dynami- 
cal systems, i.e. systems of ordinary differential equations) e.g. of IArnoldl(|1990p . 
IAmannl ( )1995p or IReitmannl p990p . The topological structure of magnetic fields 
in the vicinity of null points is described, e.g., by|Parnell, Neukirch, and Priest 



P996p and concerning the construction of ideal MHD flows, e.g. in Nickeler, Goedbloed, and Fahi^ 
(120061) . 

Our interest is to ask which topology and geometry of the macroscopic flow 
correspond to which magnetic topology and geometry in the frame of MHD. In 
contrast to the analyses mentioned in the previous paragraph we here have to 
investigate the topological and geometrical properties of both vector fields, i.e. 
that of the plasma flow and that of the magnetic field. To perform this analysis 
we solve the resistive MHD system with a linear ansatz for both vector fields, 
i.e., we perform a Taylor expansion of both vector fields in the vicinity of their 
null points, neglecting derivatives higher than the first ones. 

That implies that we allow the stagnation point (= null point of the velocity 
field) of the velocity field v to have a small offset concerning the magnetic null 
point, but neglecting here, the influence of the second derivatives of the velocity 
field. An influence would be noticeable for significant offsets, but is negligible for 
very small offsets, justifying our assumptions. We will justify this assumptions 
afterwards at the end of section [XTl estimating the influence of an offset caused 
by the gravitational force in the case of completely linear velocity and magnetic 
fields, where the second derivatives of the vector fields vanish automatically. 

Now let xq be the offset in x-direction and yo the offset in y-direction. Then 
we can express the class of linear velocity fields with the help of their first 
derivatives, i.e. Jacobians: 



( Vn 


Vi2 \ 






\V21 


V22 J 
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yo J 



v-vf=(^;- )( : - ) (1) 

in analogy to the magnetic field B 

s^55^ff"f")f!) (2) 



B21 B22 ) \ V , 

According to |Parnell, Neukirch, and Priest] (|1996p , every magnetic field of 
the form of Eq. ([2]) can be rotated such, that it is represented by a magnetic 
flux functior0 A ~ ax^ + by^, so that a standard null point of the magnetic field 
appears with a constant current density jz = jo = const around the origin, where 
the magnetic field is given by Vj4 x Cz = B. The corresponding scenario of a flare 



^This form of the magnetic flux function implies that the bisector of the separatrix angle 
should be approximately perpendicular to the photosphere, i.e. we assume the bisector of the 
separatrix angle to be symmetric with respect to the solar surface. 



SOLA: paper.tex; 4 September 2009; 8:56; p. 4 



loop, with its symmetry axis being almost perpendicular to the photosphere, is 
given in FiglH 



secondary 
plasmoids 



flare loops 




photosphere 



bisector 



Figure 1. Scenario of a flare loop with an X-point. 

The current density is 

V X B = - AA = i^Loj^ e 

such that 

Bn Bi2 
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(4) 



The different topologies of 2D vector fields are represented by two independent 

parameters, the threshold current jt, and the Eigenvalues of B, the A^'s, or j^, 
the current in z-direction 



A, 



(5) 



implying a bijcctivc relation between the variables a and b and the threshold jt 
and the actual current jz, to be precise 



-—(.3t+]z) o=— Ut-Jz 

2 2 

{a + h) jt = —{h- a) 

Mo Mo 



(6) 
(7) 
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We will use the four variables a,b,jt,jz if it is convenient, i.e. combinations of 
the four to simplify the corresponding terms in the equations. 

The Eigenvalue Xb determines the topological structure of the field and the 
geometrical shape of the field lines. For a divergence free case there are only 
three main types of such fields: 

• the case that the Eigenvalue is zero (|jt| = \jz\) corresponds to the one- 
dimensional current sheet (degenerated case). 

• the case that \jt\ < \jz\ corresponds to field lines being topological circles 
(geometrical ellipses), if jt = {a = b), then geometrical circles. All the 
cases mentioned in the last sentence are so called 0-points. 

• the case of the so called X-points, where \jt \ > \jz\- 

Only for the case in the last item a magnetic separatriccs exist. Such separating 
field lines must exist to enable magnetic reconnection. A necessary condition 
for magnetic reconnection in 2D is that the plasma flow crosses magnetic sepa- 
ratriccs, see, e.g., IPriest and ForbesI ((2000)) , [Vasyliunail pMS]) , |Cowley| ((T976)) . 



|Sonnerup| p979p . or |Sonnerup| (|1984p . The current free case is given by jz =0, 
i.e. a = —b and should be excluded, as dissipation in such a case would have 
nothing to do with the electric current. 

2.2. The basic resistive MHD equations and assumptions 

The following analysis is restricted to pure resistive dynamics, i.e. Ohmic heat- 
ing/dissipation without any other loss terms like viscosity or heat conduction. We 
choose the coordinate system in such a way, that the gravity is directed in nega- 
tive j/-dircction (the unit vector in y-direction is Cy). The basic resistive station- 
ary MHD equations in 2D are given by (following, e.g.. IGoedbloed and Poedtsl 

mm 

V • (pv) = , (8) 

p(v • V)v = jxB-Vp~ pgcy , g>0 (9) 

Eq + Y.j,By - VyBj, ^ rjjz , (10) 

7pV-v + v-Vp = (7-l)77jf, (11) 

V X B = MoJ, (12) 

V-B = 0. (13) 

Due to stationary Maxwell equations the electric field in the 2D case has to be 



constant, i.e. Ez = Eg =const see, e.g., Titov, Tassi, and Hornig (|2004p . The 
same shape of the equations occur, if we introduce a generalized non-ideal term 
on the right-hand side of Ohm's law. As the current density is constant close 
to the X-point (see Eq.[4]), the shape of an effective resistivity is that of the 
non-idealness R (where only the z-component is non-zero), writing 

V • (pv) = , (14) 
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p(v • V)v = j X B-Vp~ pgcy , g>0 (15) 

Eo + y.By - YyB, = R-e, = R,, (16) 

7pV-v + v-Vp = {-/-l)R,j,, (17) 

V X B = ^oj, (18) 

V-B = 0. (19) 

The implication is that rj = Tyoffoctivc cx i?2, as jz 7^ is constant close to the 
null point. The non-idealness resulting from two-fluid approach can be written 
for the one fluid approach, e.g., following |Braginsky| (|1965p a;. as 

1 ^ y "\ ^ "le I ^^-^^ 

ne \ / ne^ \ dt ' 



R = ^3 + —V ^ B p,] + {^)t + V • (vj + jv) ) , (20) 



where n is the particle density, e the absolute value of the electron charge, j x B 
the Hall term, Pf, the anisotropic electron pressure tensor, ^ the electron acceler- 
ation term (in our case averaged over a 'fitting' characteristic time scale, as our 

approach is quasi-stationary) , and the tensor divergence V • (^vj + jv^ describes 

the influence of electron inertia. But even this two-fluid model does not include 
all terms which may be (additionally) responsible for 'anomalous' resistivity or 
non-idcalncss, like, e.g., electron scattering by wave turbulence generated due to 
the Buneman instabihty (often considered, e.g.. lBuechner and Elkinai 120061 and 
[Karlicky and Barta[ I20(75|) . ionization-recombination terms, non-linear electron- 
ion momentum transfer terms, electron and ion viscosities, radiative losses/gains, 
and so on. As we do not specify the non-idcalncss, having introduced the effective 
resistivity, our analysis can be used for all physical models with different types of 
non-ideal terms. For a more specific view but for non-local models of reconnection 
with a generalized Ohm's law, see, e.g., [Craig and Watson] (f2003p . The term R 
is therefore only a necessary 'trigger' enabling reconnection, not a sufficient one, 



see, e.g., Schindler, Hesse, and Birn (|1988p . 

We will calculate the (almost) complete solution space of the resistive (non- 
ideal) MHD close to the null points of fiow field and magnetic field. The aim 
is to find the general correlation between the Jacobians of the plasma velocity 
and the magnetic field, and the shape of the unspecified non-idcalncss. At first 
both vector fields are treated as completely linear, i.e. unbounded fields in an 
unbounded domain, see section [3TT1 This does not exclude the possibility that in 
a bounded domain the fields at the boundary are maybe not analytic, but the 
extrapolation of the field may produce reasonable finite and regular solutions. 

In section l372l we concentrate on the linearized fields and take only first order 
terms of the spatial variables into account. A similar method to determine the 
structure of the non-ideal term and the flow and the magnetic field has been 
proposed and done, but for the case jz = and for 'global' fields, in the frame 
of a toy model in INickeler and Fahrl pOOSp . 

With this linear or linearized fields we can draw conclusions with respect to 
the other MHD quantities, like pressure, density and resistivity: The lowest order 
of the magnetic field is linear, the current density is constant, and therefore the 
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Lorentz force j x B is linear in x and y. This is the reason that all other terms 
should be also at (the) most of first order in the spatial variables. This leads to 
the conclusion that the plasma pressure p is at highest order quadratic in x and 
y, to allow for a 'linear' pressure force, i.e. we can express p by 

P^Pq+Pix + p2y + Psxy + pax^ + p^y'^ , (21) 

where po to p^ are constant coefficients. 

As the velocity close to the stagnation point can be represented as a linear 
term in x and y, the v • Vv term is also linear in x and y, or constant. The mass 
continuity equation, Eq. ([8|) implies therefore that the mass density, p, has to be 
constant in lowest order close to the stagnation point. 

Of course p could be a simple linear function close to the null point but that 
would imply that the mass density has to vanish on a surface, or rather straight 
line in 2D, depending on the gradient of the density. This null- line could be 
somewhere inside the 'linear' region. In one direction the density could then 
become negative. We assume that the density p is constant in the vicinity of the 
null points of magnetic and flow fi.eld to prevent p to have a negative or zero 
value within the domain of interest. In addition it is a reasonable assumption 
that close to a stagnation point there is a 'stagnation region' , i.e. the density 
has a maximum or minimum close to the stagnation point. 

In the case of completely linear fields it is now evident that the mass continuity 
equation requires a divergence free velocity field, if we assume all terms to have 
the same or comparable orders in x and y. 

One has also to take into consideration that at the magnetic null point Eq 
equals rjjz and if we assume that the resistivity is constant, then Eq = rjjz 
everywhere in the vicinity of the neutral point. The same holds for special cases 
of non-constant resistivity, e.g., if 77 is a function of the current density jz only, 
one can conclude that in the case of a linear field or in the vicinity of the null 
point the current density is constant and therefore also the function 77 =const, 
as 77 = ri{jz = const) =const . 

The electric field Eq is constant also outside the linear region, but the resis- 
tivity and current density are, of course, not constant there. This implies that 
the term VxBy — VyB^ equals to zero within the linear region and therefore the 
flow is field-aligned everywhere within this region in lowest order of the spatial 
variables. To get not such 'trivial' solutions, which are non-reconnective in lowest 
order, we have to regard the resistive Ohm's law as a definition equation for the 
spatially dependent resistivity, i.e. 77 = r](x, y). Therefore, to enable reconnection 
in the vicinity of the null point, the resistivity cannot be a function of the current 
density only. We designate and regard this coefficient as an 'effective' resistivity 
or short resistivity, even if this coefficient originates from collisional theory. 
As the current density is constant inside within the region of the linear field 
approximation, the (effective) resistivity is a substitute expression for a general 
non-ideal term, determining, but vice versa also determined by, the flow and 
magnetic field line structure! This relation between the velocity or flow field v, 
the magnetic field B and the resistivity rj will be analysed in this article. 

One can also recognize that Eq. (fTT|) (together with the incompressibility) 
at the stagnation point would require a vanishing resistivity in the absence of 
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additional dissipation terms. Written with all coefRcients and comparing with 
all orders of x and y we get the following system of equations, first from the 



Euler or momentum equation 

p{Vli+Vi2V2i)+2pi-2aj, = 0, (22) 

P (n'l + ■^'12 ■^21) + 2p5 - 26j, = , (23) 

P3 = 0, (24) 

p{V2i ~ V2i)ViiXf, - p{Vi2V2i+V^^)yo+P2+ pg = 0, (25) 

- pV^'-^xo - PV12V21X0 + Pi = 0, (26) 

and the equations from the energy equation, first neglecting source or heating 
terms, only including Ohmic heating, 

(-V11P1 - V2lP2)xo + {V11P2 - Vi2Pi)vo = qEq , (27) 

^21^2 + Vii (pi - 2x0^4) - 2yi2?/oP4 = 2aq{ViiXo + V12J/0) (28) 

V12P1 ~ 2V21X0P5 + Vii {2yop5 ~ P2) = 2bq {V21XQ ~ Vnyo) (29) 

21/12P4 + 2F21P5 = g(-2ayi2- 261^21), (30) 

2p4t^ii = -2agFii, (3f) 

-2p5Vn = 2bqVn, (32) 



where q := (7 — f)jz- The equations above are ordered after their physical 
meaning, the first five equations Eqs. - correspond to the first order 
momentum equations, while Eqs. (|27p - p2p represent the terms of the energy 
equation. The Eqs. (pS)) . and (P7|) are of zcroth order in x and y. The 

Eqs. (HH), (USD, (1211) and (HH) are of first order in x and y, the Eq. dSHl) is of the 
order x?/ and Eq. (^Tl) of second order in x, respectively Eqs. is of second 
order in y. 

3. Results 

In the first subsection we concentrate on linear fields, i.e. fields that are un- 
bounded in an unbounded domain. In the second subsection we take only zero 
and first order terms of x and y into account, i.e. we regard a Taylor expansion 
of maximum order one. 

3.1. Completely linear fields 

If Fii 7^ then from Eqs. (1211) and ^ it follows 

Pi = —qa (33) 
P5 = -qb, (34) 

and with two terms of the momentum equation, namely Eq. (|22p and Eq. (|23p . 
excluding the current free case jz = 0, i.e. a = —b and the case with 0-points, 
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i.e. a = b, we infer that 



p{V^\+Vi2V2i) = 2qa + 2aj, 
p{Vi\+Vi2V2i) = 2qb + 2bj^ 



Z 1 



(35) 
(36) 



thus g = (7 — l)jz = —jz should be vaUd, leading to 7 = 0. As 7 = is not 
representing a realistic equation of state, Vn must be zero, taking terms of every 
order into acount. 

This is unaffected by the neglecting other source terms, not mentioned here, 
as the interesting heat conduction term is of the order zero in x and y, and 
therefore affecting the Eq. (|27|) only. 

Now we use the assumption that Vu = 0. Then we rewrite the MHD equa- 
tions: 



pVi2V2i+2p4-2a], = 0, 

pVi2V21+2p5~2bu = 0, 

-PV12V21 j/o +P2 + pg = , 

-pVi2V2ia;o +Pi = 0, 

-Vi2yoPl - V2lXt^P2 = qEq , 

V21P2 - 2Vi2j/oP4 = 2qVi2ayo , 

V12P1 - 2V'2ia;oP5 = 2qV2ibxo , 

Vi2{p4 + qa) + V2i{p5+qb) = 0. (37) 



With Eqs. ([To]) and (fTTjl it can clearly be recognized that at the stagnation 
point the resistivity must be zero. As the electric field Eq is constant and Eq ~ 
rjjz at the magnetic null point and the stagnation point, we conclude that the 
electric field is zero (everywhere) . The effect of a vanishing compressibility is in 
resistive MHD also connected with a vanishing electric field. 

In addition, as the x-offset xq in the fourth equation (a zeroth order term of the 
force) of the system ([37]) couples an additional pressure gradient in a;-direction to 
this offset only, this offset can be set to zero without loss of generality. Therefore, 
from xo = and Eq = we infer pi = 0, see the fourth and fifth equation of the 
system (j37p (part of the energy/entropy equation). 

If we would further allow yo — in the case of non-vanishing gravity, then 
there would exist only static solutions either (as in this case V21 must be zero 
it implies that V12 = 0), or the other possible solution branch would need an 
adiabatic exponent 7 = 0. 

For the non-degenerated case {yo 7^ 0) the general solution for the system 
([37]) has two branches: 



P2 = 



pg _ p^gVi2 

2 ib-fj. 
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P5 = U{2~j)bj. + P^±lVD 



yo 



P9 



= ^^-TT , (OOj 

2pVl2 

where the discriminant D must fulfill 

D = 8apVi|7J. + - pVi^ ' > (39) 

to guarantee that V21 is real. This leads to 

D = {pV,\f + (8a - 46) ^j,pV,\ + Ab^j^jj > (40) 
^ i?2:=[(4a-2&)7j,]'-46Vj'>0, (41) 

where D2 is the discriminant of the quadratic equation Eq. (|40[) . The inequality 
Eq. ((1T|) guarantees a real value of pT^fl • That leads together with the convention 
Jt > 

(a - 6) a < ^ > -Jt . (42) 

This implies that either the actual current density jz > 0, then the field lines 
are hyperbolas or topological circles, or jz < 0, then only hyperbolic field lines 
or (anti-)parallel are possible. 

Let us now estimate the influence of the offset in y-direction on the solution, 
e.g., for the situation in the corona. We can express the offset in y-direction by 
using the equation for the offset yo in Eq. ((38|) 

assuming that 7 = 5/3. 

We demand that the configuration should really be an X-point, but not a 
one-dimensional sheet, therefore we assume that the relation between the two 
characteristic currents is given by a value not much larger than jz/jt ~ 0.999 
(correponds to an opening angle of the smaller separatrix angle of only about 2 
degrees). With the electron charge e, the number density n and the drift velocity 
Vd and jz ~ nevd, it follows for yo 

yo = 1200 ^ = 1.2x103^^ (44) 

The value |yo| must now be compared with the typical lengthscale and values of 
a magnetic field close to a null point, i.e. the coronal magnetic field is about Be ~ 
10~^T (the index 'c' stands for 'coronal') and the scale on which the magnetic 
field varies linearly should be larger than a Debye length (for coronal parameters 
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in the range of mm to cm) but much smaUer than the active region loops 
{Ic ~ 10^'^m) to get an upper hmit for |yo|- We therefore use the approximate 
current sheet thickness d. To estimate the current sheet thickness of a finite 
current sheet, we assume again that the current density jz = nevd, where is 
the drift velocity, which should not exceed the thermal velocity (of the electrons), 
i.e. Vd < y/ikBT/m^. Let ks = 1-38 x 10~^^JK~^ the Boltzmann constant, 
£o = 8.85 X lO^^^A^s'^kg-^m-^ the dielectric constant, nip = 1.67 x lO^^^kg the 
mass of the proton, rrie ~ 9.1 x lO^^^kg the mass of the electron, n = 1.6x 10^'^m'^ 
the number density, and Td = lO^K (as the thermal temperature is about T — 
2 X lO^K, we consider an effective temperature Td, thus that the corresponding 
drift velocity Vd is definitely below the critical drift velocity given by temperature 
Vj^ ikBT/nie to guarantee a stable stationary current, instead of an kinetic 



instabihty see, e.g., Papadopoulos I1977P the approximate temperature (for the 



coronal parameters see, e.g.. IStixl I2002p . We choose a 'small' current density to 
get the upper limit of y^. The offset is then given by 

, , 1.2 X 10^TO„5 1.2 X 10^ m„o ^ _,„_f; 

2/0 = : ^-TT ~ 3 X 10 ^n . 45 



The lower bound of the current sheet thickness is approximately given by 

^^^=«33m. (46) 

Mojz fioney/3kBTd/me 

The fraction of both values should be yo/d land is then determined by 

Wn o 771„(7 1.2 X lO^TO^O o 

— = 1.2 X 10^ — ^ = w 6 X 10"* < 1 . (47) 

d evdBc e^3kBTd/meB, 

Thus for typical coronal parameters the offset in j/-direction has a value of about 
mm or even some orders of magnitude less, in contrast to the value of the typical 
lengthscale of the linear X-point region which is at least of the order of one meter 
or some tenth of meters or even more. This is an argument or rough justification 
to neglect the term of the interaction between the shift and the second derivative 
of V. 

Even if we assume that the drift velocity v^ is, let us say six orders, smaller 
than the thermal velocity, the fraction |yo|/'^ would be of the order of some 
percent. Therefore, even in this extreme regime of the drift velocity, the stagna- 
tion point is located within the current sheet, assuming that the sheet is really 
two-dimensional (non-singular) and its extensions in both coordinate directions 
arc about of the same order. The complete mentioned discussion is valid exactly 
only in the case of unbounded fields, but it shows that in general the offset due 
to gravity can be neglected. 



Ad = (neVCeofcsT))^'^^ 10 
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3.1.1. Discussion of the discriminant D 



The discriminant D2 in Eq. (j41[) . being larger than zero is only a necessary 



criterion. That ^1^2 > has also to be guaranteed. 

To fuUfil the relation Z? > we define two values, namely 



ei := pVkcrUi = Y (3jt + 3.) 13z + 2fiol\j. \ \l (48) 



' jt (jt 4 


-jz) 


2 







to represent the relation Eq. ([40|) 

(pVil _ (^y^2 _ ,2) > (50) 

^ C(ase)I pV^2 > ei A p^il > £3 (51) 
is valid or 

=> C(ase)II pV^2 < ei A pV\| < £3 , (52) 

which implies for CI pVi2 > max(ei, €2) = ei and for CII pVi2 < niin(ei, €2) = £2. 
The values £1,2 are explicit functions of the current densities jt and jz - This gives 
us restrictions for V12: for > the lower bound for V12 is given by the left 
relation of Eq. (|5ip otherwise the lower bound is zero, writing 



pV,l > ^{ijt+J.)lJ. + 2pol\j.\ f'^''^''\ (53) 
such that the non-admissible region is then given by 

V12 ^ [+V12 -V12 .critl ] , (54) 

and V12, critl is given by 



V12 



jtjjt+jz) 

2 



(55) 



For jz < all values of V12 are allowed. There are no restrictions for V12, as for 
jz = the critical value Vi2,criti = and for jz < all values of Vi2,criti are 
imaginary. 

For CII it is clear that the criterion is given by 



pV,l < ^{3jt+jz)lJz-2pol\jz\ f'^''^''^ . (56) 

For jz < no value of V12 can satisfy the relation Eq. ((56|) . Only for jz > 
solutions exist in the region beyond the curve Fig. [2] and the j^-axis. 
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Thus we can summarize: For jz > the ahowed values for pVi2 are given by 
the relation 

£2 > pV^2 > V pV^2 > ei ■ (57) 



3.1.2. Solutions with shear flows 



In the case of a vanishing V12 or V21 the equations again lead to the statement 
that either V12 — V21 = in both cases or that again 7 = 0. This implies 
that the allowed flows are of higher orders with respect to the spatial variables, 
like indicated in Priest and Cowle"y| (|1975p , or need special, maybe unphysical 
thcrmodynamical constraints. Also the case of general shear flows leads to 7 = 0. 
Such solutions will not be discussed here. 



3.1.3. The resistivity 

The resistivity is given by the z-component of v x _B and results in the quadric 
77 = - {2aVi2yox - {2aVx2 + 2&F21) xy) . (58) 

Neglecting the term with yo, the above equation shows that the isocontours 
of the resistivity have four branches. In two quadrants 77 has a positive sign and 
the other two quadrants represent an effective negative resistivity. 

This example of unbounded fields shows an interesting behaviour because 
the solutions are exact analytical ones and they exhibit the fact that there is a 
special relation between the derivatives of v and of B connected about algebraic 
varieties which generate the solutions. We stop these investigations here due to 
the problem of negative resistivity. To interpret this negative effective resistivity 
(or negative non-idealness) and search for a reasonable physical interpretation 
or correponding physical mechanism needs further investigation. 



3.2. Complete linearization of the resistive system concerning x and y 

Here Vn 7^ in general, as the equation which contradicts the thcrmodynamical 
constraints, the last two equations of Eqs. (|5T|) and ((5^ . being part of the energy 
equation, are neglected due to their higher order in x, y. 

In this case we only need to regard equations, being of first order in x and y, 
namely 



P {V,\ + V12V21) + 2p4 - 2ajz = , (59) 

P {V^\ + V12V21) + 2p5 - 2bjz = , (60) 

P3 = 0, (61) 

~p{Vi2V2i + V^^)yo+P2 + pg = 0, (62) 

-pV^j^xQ- PV12V21X0+P1 = 0, (63) 
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Figure 2. Not all values for the 'free' velocity parameter V12 are allowed. There exist some 
forbidden regions for the velocity parameter. 



(- Viipi - 1^21^2)2:0 + (V11P2 - Vi2Pi)yo = , 



(64) 
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V21P2 + Vii {pi - 2.T0P4) - 2V12J/0P4 = 2aq{ViiXo + Vuyo) , (65) 
V12P1 - 2V^2ia::oP5 + Vn (2yoP5 - P2) = "^.hq (V21.T0 - Vnj/o) , (66) 



which contain only terms of first order. This is a nonhnear system consisting of 
eight, but. of course, effectively seven equations for nine unknowns. 

3.3. Solutions without gravity and shift 

In this case the stagnation point is identical with the magnetic null point. The 
general solutions is now represented by 

P4 = (V'n + V12F21) + aj, 

P5 = -^{V^i+Vi2V2i) +b3,, (67) 

with free parameters Vii and V12 (and, of course, the 'magnetic parameters' the 
threshold current jt and the actual current jz). The parameter V21 is not free, 
as it is determined by the fact that the quadric (surface), representing the 
resistivity, must be an elliptic paraboloid. We have to introduce another pa- 
rameter s to express V21 as function of V12 and s. The elliptic paraboloid is the 
only quadric that allows a positive resistivity with one zero (the vertex of the 
paraboloid) at the magnetic null point. All other quadrics can be excluded, 
with the exception of the degenerated case of a cylindrical paraboloid, see, 
e.g., Bronstein and Semendjajew p987p or IBartschI (|1984p . The cylindrical 



paraboloids are limiting cases with s = ±1 as is shown as in the example in Fig.[3] 
(here the specific value is s = ^1)- The general solution is then parameterized 
by Vii, V12, s 



1^21 = ^^^^^Vi2 + 'isVii , (68) 

jt - jz jt - jz 

with the restriction s e [—1,1]. Multiplication of Eg. ([55]) with V12, inserting 
this expression into the expression for the Eigenvalue Ay = + V12V21, and 
completing the square gives 

\l = V^^+Vi2V21 



= 1^11+5^12^^^^-^ 1 +(l-s2)V,|^^>0. (69) 

Y Jt- Jz J Jt- Jz 

If \jz\ < jt and \s\ < 1 as assumed, the magnetic field is of hyperboHc type and 
it is guaranteed that the stagnation point is also of hyperbolic type, i.e. has two 
purely real Eigenvalues, namely one positive value and a negative counterpart. 

Although the solutions shown here are no reconnection solutions, as the elec- 
tric field has to be zero, because only resistive dissipation is included in our 
investigation, they represent almost reconnection solutions(=reconnection-like 
solutions) where the necessary condition for reconnection is fulfilled. 
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Figure 3. Flow for the parameters Vii = V12 = 1 and s = (left column) and s = —1 (right 
column), crossing all four magnetic separatrix branches (top panels), shape of the positive 
resistivity (middle panels) , and the elliptic shaped contours of the resistivity (bottom panels) . 



Let US investigate now the geometry, respectively the slopes of the separatrices 
of the flow and the magnetic field. Defining Jz := jz/jt, the magnetic separatrix 
is given by A = and A ^ ^ [{jt — jz)y'^ — (jt + jz) x^)] (we restrict this to 
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\Jz\ < 1) and can therefore be expressed by 



V^-JzV ± V'TTX.T = ^ y = ± ^^^'^' x := ±KmX , (70) 

V 1 ^ Jz 

where ±Km is the slope of the both magnetic separatrix Hnes or to say the 
both asymptotes. As V21 is given by the general solution. Eq. ([68|) . by defining 
a stream function tp via 



dip , dip 
— and V =-—, 
ay ox 



— and Yy = ~—, (71) 



we can integrate the above equations and get for the stream function 

^^^V'~^x' + V,,xy. (72) 

The fluid separatrix is here given by -0 = 0, geometrically this arc asymtotical 
lines (asymptotes). 

We will briefly discuss the problem that the magnetic separatrix is partially 
identical with the hydrodynamic separatrix. In this case the plasma flow can 
only take place across one part of the separatrix, or both separatrix lines are 
identical and no rcconncction can take place. 

One can clearly recognize that in the case Vn = and thus {V21IV12Y = 
(1 + J^)/(l — Jz) both corresponding asymtotical branches(=separatrix lines) 
have the same slope. Thus the hydrodynamical separatrix and the magnetic 
separatrix are identical, the plasma cannot cross the magnetic separatrix and 
therefore no reconnection can take place. 

Writing y{x) = Kyx (if Vii,Vi2 7^ 0), where Ky is the slope of the hydro- 
dynamic separatrix, and inserting this into -^y^ — ^-^x'^ + Viixy = G with the 
parametric expression for V21 in Eq. (|68p. we get the slopes Kviy2 of the two 
asymptotes / separatrix lines 



Vii ^ V^^+Vi2V21 



V^2 



Vrel ± W + 2sVrel \ _ j + J— ^ , (73) 



l + J. 

l-Jz 1 - J; 

where Vrei = Vn/Vu- For V12 ~ ^ V\\ the asymptotes are given by x = and 
y = 3-^(1 + Jz)/(1 — Jz)x, and for V\2 — V21 = we get x = and y = as 
separatrices. 
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Let Vii = Vl2(= 1), then 



Kvia = -l± 



I -J. 



= -l± 1 



l-Jz 



(74) 



There is one corresponding slope of the flow separatrix with respect to one of 
the magnetic separatrix hnes, if s = ±1. 

Therefore we can formulate the following theorem (in analogy to known anti- 
reconnection theorems) , which is now restricted to almost reconnective solutions 
with vanishingly small electric fields and resistive disspation only, i.e. vanishing 
non-resistive dissipation mechanisms: 



Theorem If the flow close to the null point is in good approximation incom- 
pressible, the electric field is negligible and other kinds of dissipation mechanisms 
than the resistive dissipation are also negligible, then the plasma flow cannot cross 
the magnetic separatrices if for the both slopes of the both magnetic separatrix 
lines K„i = Kmi = —Km2 with Km = ■\/(l + Jz)/(1 — Jz) and the following is 
valid: 



I. a) Km = Kvi A h) - Km = Kv2 
or 

II. a) Km ^ Kv2 A b) - Km = Kvi (75) 



// only I. a) or Lb) or II. a) or lib) is valid then only one magnetic separatrix 
can be crossed and as the other stream lines converge to the second magnetic 
separatrix line without crossings, one could call this reconnective annihilation, 
see, e.g., Priest & Forbes (2000). The necessary condition for a complete non- 
crossing is Vii = for Vi2 ^ and for only partly crossing flows s = ±1. 

For almost all values of s, here s = 0, the flow in Fig[3] (left column) crosses 
all four separatrix branches (all two separatrix lines), the resistivity is positive 
(only zero at the null point, middle panel), and their isocontours are ellipses 
(bottom panel). For the special parameter s = — 1 it can be seen in Fig. [3] 
(right column) that the flow, like for the aforementioned magnetic reconnective 
annihilation solutions, crosses only the separatrix line with the positive slope, 
while it converges to the other magnetic separatrix, being almost field aligned 
(top panel) , the resistivity being also postive (middle panel) , and the isocontours 
of the resistivity are straight lines (bottom panel). 

We now focus on the effective resistivity and therefore on the nature of ohmic 
heating/disspation. The resistivity is given by 



Mo Vii 



[jt + jz) - 2s^if - jlxy + {it - jz) y" 



(76) 



To ensure that the dissipation is positive we prove that the term {jt + jz) ~ 
1s\J j1 — j^xy + {jt — jz) y^ is larger than zero (or zero) for all s and x, y. Let 
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Ki^2 = jt i jz and Ki, K2 > as required for an X-point. Then 



<Kix^ -2s^/K^xy + K2y^. □ (77) 

This imphes that the sign of the resistivity rj depends on the sign of -yJ- . As the 
dissipation should be positive, i.e. 77 > 0, thus Vn and jz must have the same 
sign. 

For s = ±1 it can be recognized from Eq. ([75)1 and from the bottom panel 
Fig. [3] (for the case s = —1) that the resistivity contour lines are straight lines 
given by 77 = ^o^ii/(2j^) (Vjt + jzX ± VJt^^y)'^, and therefore the resistivity 
has the structure of a one-dimensional sheet. The case s = ±1 therefore reflects 
the ID character of the non-idealness, which can also be found in the papers of 
[Craig and Henton| (fTMS]) . 

4. Conclusions 

We analyse the solution space of the MHD system with a generalized non- 
idealness in Ohm's law close to an X-type magnetic null point. This procedure 
is done to get non-ideal/resistive but formally non-reconnective solutions that 
can be regarded as reconnection solutions with a vanishing electric field. Such 
solutions show the parameters and their relation amongst each others, to get non- 
ideal non-reconnective, i.e. either almost reconnectivc annihilation or completely 
almost reconnectivc (= reconnection-like) solutions. 

We assume that the flow close to the magnetic null point is governed strictly 
by the pure non- ideal/resistive system of MHD equations, and that the stagna- 
tion point is close to the null point (or that both singular points in the plane 
are identical). The focus is on the region close to the magnetic null point to 
investigate the possible type of topology of the corresponding stagnation point 
(flow). To use a closure for the equations and to analyse the resistive MHD 
system we use the corresponding energy equation. 

We further assume that the density is constant close to the stagnation point. 
This assumption is not a necessary, but a plausible condition concerning the 
mass distribution in the vicinity of stagnation points. This leads to a vanishing 
divergence of the plasma velocity and as a logical consequence to a vanishing 
resistivity close to the stagnation point. Inserting this into resistive Ohm's law 
implies that the electric field has to vanish at the stagnation point. We conclude 
that the electric field has to be zero, and as a consequence of the 2D stationary 
approach the electric field has to be zero everywhere. 

This case of vanishing divergence of the velocity field and thus of a vanishing 
electric field can be regarded as a physical approximation or limit. Obviously a 
classical reconnection process requires a non- vanishing divergence of the velocity 
field to produce a non-zero reconnection rate, determined by the constant electric 
field Eq, see, e.g.. lPriest and ForbesI (pOOO]) . 

From resistive Ohm's law one can clearly recognize that close to the magnetic 
null point the electric field is identical to the 'convcctive' electric field generated 
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by the product from resistivity and current density. As the electric field is zero, 
the resistivity inside the finite current sheet must be zero at the magnetic neutral 
point. The effect is that the resistivity cannot be constant close to the magnetic 
neutral point. Therefore the resistive Ohm's law is an equation determining 
the topological type of the magnetic field and the velocity field and connect- 
ing/relating it to the shape of the resistivity, which is a quadratic function of the 
spatial coordinates. We call this spatially varying resistivity effective resistivity 
or for short resistivity as usual. As the current density is constant close to the 
magnetic null point, the right hand side of Ohm's law, Eq. (|10p could be regarded 
as a generalized non-ideal term. 

A necessary condition for reconnection is that the plasma flow can cross 
magnetic separatrices. Wc found that crossing of magnetic separatriccs requires 
also an X-type stagnation point flow if the effective resistivity and thus the 
dissipation should be positive. 

The sufficient criterion for reconnection is the crossing of the separatrix and 
the non- vanishing electric field. This reconnection process can take place for/in 
the incompressible case if there is not only a resistive, but an additional non- 
resistive dissipation (term) that occurs only in the energy equation. Therefore 
one possibility could be heat conduction. Another possibility is the viscous case 
with constant viscosity. Here, the Laplacian of the velocity field vanishes in the 
first order momentum equation of the ions (Navier-Stokes) as the velocity field is 
linear. These cases have to be investigated in more detail but have been discussed 
under the assumption of constant resistivity, e.g., by IPriest and Forbesi pOOOp . 

The analyses started by us here will be extended to configurations that have 
density gradients close to the stagnation point and therefore allow for a non-zero 
electric field. The extension to energy equations allowing for a deviation from 
the classical energy equation of resistive MHD will provide us with a relation 
between the deviation from pure resistive dissipation and deviation from the 
classical (X-type) stagnation point flow. 
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